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Abstract 

We consider the Kepler problem on surfaces of revolution that are homeomorphic 
to 5"^ and have constant Gaussian curvature. We show that the system is maximally 
superintegrable, finding constants of motion that generalize the Runge-Lentz vector. 
Then, using such first integrals, we determine the class of surfaces that lead to block- 
regularizable collision singularities. In particular we show that the singularities are 
always regularizable if the surfaces are spherical orbifolds of revolution with constant 
curvature. 
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1 Introduction 

The problem of describing the motion of a particle on surfaces of constant 
curvature, under the influence of the analogue of the gravitational potential, 
is an interesting problem that dates back to the 19th century. A short his- 
tory of the problem is presented in [10]. Many of the classical results have 
been long forgotten. However, since then, interest in the problem has been 
revived, at least in part, because of cosmological models as the one based on 
the Robertson- Walker metric. This model describes a homogeneous, isotropic 
expanding or contracting universe, and its spatial slices are, depending on the 
curvature, three-spheres S^, copies of the Euclidean space M^, or copies of the 
hyperbolic space H^. In recent years many authors have studied the classical 
Kepler problem and the quantum analogue (the hydrogen atom), rediscovering 
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the old results and introducing new ones (the interested reader can consult [2] 
for some interesting results and an extensive bibliography) 
In this paper we study the Kepler problem on surfaces of revolution of constant 
Gaussian curvature with certain type of metric singularities. This problem is 
more general than most of the previous work that concentrates on constant 
curvature spaces with no singularities (i.e. the sphere 5*^ and the hyperbolic 
plane H^, or in higher dimensions the three sphere 5*^ and the hyperbolic 
space H^). 

As in the standard Kepler problem the potential is singular and this introduces 
singularities in the equations and in the solutions. In our case we will consider 
spherical surfaces of revolution that have constant curvature, i.e. surfaces of 
revolution that are homeomorphic to the two-sphere S'^ and that can have 
metric singularities at the "north" and "south" poles. 

Levi-Civita has given us a method for extending orbits of the Kepler problem 
through double collisions [7] . He eliminated the singularities in the vector field 
by transforming the equations to ones without singularities. The extension 
through double collisions then is given by the transformed equations and is 
automatically a smooth function of initial data. In [6J this method is applied 
to the case of the Kepler problem on spaces of constant curvature (with no 
metric singularities). 

Another method of regularization is the topological regularization due to Eas- 
ton [5], that is based on the general theory developed by Conley and Easton 
and it is usually called regularization by surgery, or (following McGehee 
[9J) block regularization. This method is used in this paper to regularize the 
singularities at the poles. Roughly, the idea is to excise a neighborhood of the 
singularity (more precisely an isolating block) from the manifold on which the 
vector field is defined and then to identify appropriate points on the boundary 
of the region. This is done constructing a map across the block that identifies 
the point where a solution enters the block with the point where a solutions 
exits the block. This map is a diffeomorphism. However there are solutions 
that, once they enter the isolating block, never leave it. If the summentioned 
map can be extended (in a differentiable way) to consider such solutions then 
one says that the singularity is block regularizable. 

We characterize, as in [10] , the spherical surfaces of revolution using a param- 
eter (3. In [TU] we showed that whenever /? is a rational number the Kepler 
system has the Bertrand property (i.e. all the bounded non-singular orbits 
are closed) Here we find that, while the south pole singularity is always block 
regularizable, only few values of P produce a singularity at the north pole 
that can be regularized according to Easton. Furthermore, it turns out that 
the north pole singularity is always block-regularizable in the case of closed 
(compact, without boundary) surfaces of revolution (with constant curvature) 
that are Riemannian 2-orbifolds. 

Loosely speaking, a 2-orbifold is modeled locally by convex Riemannian sur- 
faces modulo finite groups of isometrics acting with possible fixed points. This 
means that a neighborhood of each point p of such orbifold is isometric to 
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a Riemannian quotient Up/Tp where Up is a convex Riemannian surface dif- 
feomorphic to M^, and Tp is a finite group of isometrics acting effectively on 
Up. Every Riemannian surface is trivially an orbifold, with each Tp being the 
trivial group. The reader interested in more background on orbifolds can con- 
sult [llj. For the purpose of this paper, however, we only need to apply a 
simple explicit criterion to determine whether a closed surface of revolution is 
a 2-orbifold (see section 10 and JT\). 

The paper is organized as follows. In Section 2 we introduce the the general- 
ized Kepler potential and the equations of motion. In Section 3 we introduce 
Gaussian curvature and some properties of constant curvature surfaces. In the 
following section we find two additional integrals of motion, besides the Hamil- 
tonian and the angular momentum. In Section 5 we find the equations of the 
trajectory. In Section 6 we use a transformation to rewrite the equations in 
a more convenient way and we define the collision manifold. In the following 
section we study the flow on the collision manifold. In Section 8 find that for 
only few values of (3 there is a trivializable isolating block about collisons. In 
Section 9 we determine for which surfaces the singularities are regularizable. 
In the last section we show that the singularities are always block regularizable 
if the surface under consideration is an orbifold of revolution. 



2 Equations of Motion 

Let I be an interval of real numbers then we say that 7 : / — > is a regular 
plane curve if 7 is and 7'(x) 7^ for any x G /. 

Definition 1. Let 7 : / — be a simple (no self intersections) regular plane 
curve 7(r) = {f{r),0,g{r)) on the x^-plane where / and g are smooth curves 
on the interval / = [risi,rs], with /(r) > and /(r) = if and only if r = tat or 
r = r^. A spherical surface of revolution is a surface isometrically embedded 
in that admits a parametrization x : / x M ^ S* of the form 

x(r, 0) = (/(r) cos 6, f{r) sin 6, g{r)) (1) 

That is, S is the surface of revolution obtained by rotating 7 about the 2;— axis. 
The curve 7 will be called the profile curve. 

Note that a spherical surface of revolution is homeomorphic to 5*^ and that by 
definition the sets x(rAr,6') and x(r5,^) reduce to single points, i.e. the north 
and the south poles of S. Metric singularities can only occur at the north and 
south poles, S is smooth everywhere else (see figured]). 

Throughout this paper all surfaces of revolution will be assumed to be as in 
Definition 1 and the profile curve 7 is assumed to be unit speed, i.e. (^)^ + 

(d£\2 _ 1 

For a surface of revolution S, a simple computation gives the coefficients of the 
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Fig. 1. A spherical surface of revolution with positive constant curvature 
first fundamental form, or metric tensor (subscripts denote partial derivatives): 

£; = x,.x,= (^^j =1, F = x,.x, = G = X, • X, = /(rf , 

so that the metric (away from any singular point) is 

ds^ dr^ + 2F drde + G dO^ = dr^ + f{rfde\ (2) 

Note that the parametrization is orthogonal {F = 0) and that Eg = Gq = 0. 
Surfaces given by parametrizations with these properties are said to be r- 
Clairaut. The Lagrangian function of a particle of mass m moving on the 
surface takes the form 

where V{r) is the generalized gravitational potential, that is 

V(r) = 7e(r) (3) 

where 7 is a positive constant and 0(r) is an antiderivativc of l//(r)^. The 
generalized gravitational potential is a solution of the Laplace-Beltrami equa- 
tion 



The Hamiltonian is 



where pg — mf{rY9. 
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The equations of motion are: 



dH 



Pr_ 

m 



= ^ = Pfl 

dps mf(rY 
^ _dH. = Pef'i''') 



dV 
dr 



(5) 



Clearly H and pg are constant of motions, they are in involution and the 
problem is integrable by the Liouville- Arnold theorem. 



3 Curvature 

It is well known that the (Gaussian) curvature of the metric is given by 
the equation 

fir) 
fir) 

Throughout this paper all the metrics will be assumed to be of constant cur- 
vature K. The requirement of constant curvature gives us a linear differential 
equation to solve 

/" = -Kf. 

The solutions to this differential equation are of the form 

/(r) = Aoe'"^' + Boe-'^'' 

ifK^O and 

/(r) = Cr + D 

UK = 0. 

Since we are interested only in spherical surfaces of revolution we can restrict 
our attentionto the case of positive curvature, i.e. i^' > 0. In this case the 
solutions take the form 

f{r) = Asm{y/Kr)+Bcos{VKr) or f (r) = L sm{\/Kr + a) 

where A = L cos a and B = L sin a with L > 0. 
Now we can prove few useful facts and formulas. 
Proposition 1. The equation 

- //" + iff = (6) 

is verified if and only if the metric has constant Gaussian curvature K and 
either /(r) = Aqc'^'' + Boe-'^'' with AoBo = or /(r) =Cr + D with 

C = ±f3. Moreover if the curvature K is positive, then A"^ + B'^ = (3"^ / K. 
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Proof. Note that 



\ P J f 
If -//" + iff = P'^ then from Eq. (CD it follows that 

fif?~P'^ 



f 



-K 



for some constant K. Consequently, since —ff" + if'Y = /"// = —K and 
the curvature is constant. 

On the other hand assume that /" = -Kf. Then, if /sT ^ 0, / = Aqc'^'' + 
^^g-iv^r_ Plugging this into -//" + {ff = p'^ we find the condition AqBq = 
jj^. If i^' = then / = Cr + D. Plugging into the equation we find = 

□ 

Proposition 2. The function / satisfies the equation 

fir) 



fir) 



-(3'e{r) 



for some antiderivative 0(r) of l//(r)^, if and only if it satisfies the nonlinear 
differential equation 

-//" + if)' = 



Proof. 

(f\ _ f"f-ifr _ 
[fj p p' 

which implies (|8]) for some 0(r). □ 

Proposition 3. If / satisfies the equation — //" + (/')^ = and 6(r) is 
defined as above, then 

Proof. From Proposition [2] and Proposition [T] we obtain 

if)' P' + fP 



Since, by Proposition [T], the Gaussian curvature is constant we can substitute 
/" = —Kf in the last equation. The result follows. □ 
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4 Additional Integrals 



In Section 2 we found that the system under consideration has two integrals of 
motion: the Hamiltonian and the angular momentum. In this section we con- 
sider the Kepler problem on surfaces of revolution that have constant Gaussian 
curvature and we find two additional integrals of motion (quadratic in the mo- 
menta). Three of the first integrals are independent and thus the system is 
maximally superintegrable. 

We now look for first integrals that are quadratic in the momenta. The most 
general form of an invariant, quadratic in the momenta is 

I = apl + 2bprpe + cp] + $(r, 9) 
where a, 6, c and $ are function of r and 6. Given the form of the constant of 
motion , it is straightforward to find the conditions that must be satisfied to 
grant its existence. The Poisson bracket of / with the Hamiltonian is 



'dIdH dl dH\ f dl dH dl dH" 



\ dr dpr dpr dr J \d9 dpe dpg dO J 

3^(oh^ «e \ 2 ^( 2ar 2he\ , (2bf' eg' 

--arP, + ^2br + j p^pe +[Cr + J^ + J^j PrPe + y-jr + J2^ 

+ - 2aVr)Pr + \-2hVr 




where we used the fact that V is a function of r only and thus Ve = 0. 
For J to be a constant of motion the Poisson bracket {/, H} must vanish for 
every value of the momenta. The vanishing of the Poisson bracket implies the 
set of equations 

Ir = 



2^ + £| = 

f ^ (9) 
26^ + M = 



and 



<l>^ - 2aVr = 0, -2bVr + ^ = (10) 

/ 

We now consider the motion under the generalized gravitational potential and 
we look for additional integrals of motions. Thus let V{r) = 7G(r) then the 
equations (fTOl) take the form 

= 2^, $e = 27&. (11) 
The compatibility condition for the last two equations leads to 
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or ae = p{r)hr- On the other hand from ([9]) we have ag = — 26.f./^(r) and thus 
br = and b = b{6). Moreover by (Q we also have a = a{6) and, since 6^ = 
we have that ag = 0. Consequently a = constant. 
Integrating the last equation of system we obtain 

c{r,9) = ^-2bg^ + m 
f 7 

where F is an arbitrary function of 6. Differentiating with respect to 6 yields 

V(r) 

eg = -2beg^^ + Tg. 

7 

On the other hand the second equation of system iQ implies 

c. = -2-. 

Therefore we obtain 

V f 
-2bge- + Tg = -26^. 
7 / 

Let V = 79(r) and let — //" + (/')^ = then by Proposition [21 we obtain 
or if we assume F = constant we obtain. 

bgg + P% = 0. 

The last equation is the equation of an harmonic oscillator and its general 
solution is b{6) = Eism{i36) — -E2 cos(/56'), where Ei and E2 are arbitrary 
constants. Thus we find a constant of motion / of the following form 

I = 2aH + Vpl + 2^1/1 - 2E2I2 

where H is the Hamiltonian and 



h = sm{p9)prpg -Pcosipe)e{r)pl - ^ cos(/3^), 



l2 = - cos{/3e)prPg - l3sm{(3e)Qir)pl - ^ sm{(3e) 



(12) 



Note that 2aH + Tp^ is a constant of motion and, since Ei and E2 are arbi- 
trary constants, Ji and I2 are integrals of motions that generalize the Laplace- 
Runge-Lenz vector. 

The four conserved quantities H,pg, Ii and I2 are related by the equation 

l! + l! = 2pjH-^ + l-^ (13) 
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giving three independent constants of motion. Equation (fT3|) can be easily 
derived from ^ 

using Proposition [3l i.e. using that 6^ = -^ji — ^- Substituting the expression 
for in the previous equation we obtain 



5 Equation of the Trajectory 



The shape and orientation of the orbits can be determined using the general- 
ized Laplace-Runge-Lenz integrals: 

Ji cos(/3^) + h sin(/?^) = -/5e(r)pe' - ^ 



= ^ I 1 + -(cos(/3^)/i + sin(/?^)/2 



and thus 

n V 7 
or 

p= -(l + ecos(/3(0-0o)) (14) 

where p = -Q{r), p = (^-^^ , (hP/l) = ecos(/36'o) and (hP/l) = 
e sin(/?^o)- 



6 A Geometric Description of The Flow 



From now on we consider only the Kepler problem on spherical surface of 
revolution with positive constant curvature. The sphere is the trivial example 
of spherical surface of revolution with positive constant Gaussian curvature. 
All the other surfaces of revolution of this kind have metric singularities at 
the north and south pole. An example is depicted in figure [H 
We now present a description of the orbit structure of the system, with special 
emphasis on the orbits near collision. The coordinates used here are a gener- 
alizations of those of R. McGehee used by several authors to study collisions 
in Newtonian gravitational systems [S]. If the Gaussian curvature is constant 
and positive, by Proposition [21 since /(r) = Lsin(v^r + a) and = L^K, 
we find 

e(r) = = cotfv^r + a). 

Clearly /(r) = when sin(V^r + a) = 0, i.e. when ^J~Kr + a = or when 
^/Kr + a = vr. Let r^r = —a/\/li and r5 = (vr — a)/\/K. Then /'(r) = 
L\/Fcos(v^r + a), thus /'(r^v) = v^L > and f'{rs) = -VKL < 0. 



9 



Consequently V{r) has an attractive singularity at the north pole r = and 
an equal repulsive singularity at the south pole r = r^. Since V : (r^r, r^) M 
is real analytic, standard results of differential equation theory guarantee, for 
any initial data {r{0),9{0),pr{0),pg{0)), the existence and uniqueness of an 
analytic solution defined on a maximal interval [0,t*), where < t* < oo. If 
t* < oo, we say the solution is singular. In general there are different kinds of 
singularities of the solutions. However in this problem we have only one kind: 
the singular solutions are such that r{t) — »• rjy a.s t —>■ t*. In this case we say 
that the solution experience a collision. Thus the singularity of the potential 
at the north pole r = rjy induces singularities in the solutions and corresponds 
to a collision (on the other hand the singularity of the potential at the south 
pole r = rs does not induce singularities in the solutions). 
To study the flow near collisions, i.e. near r^v consider the following transfor- 
mation of coordinates 



Pr 



U - 2^ 



In these coordinates, taking m = 1, the original system becomes. 



dr 
dt 

<fo _ (,,2 f,(„\ _ 2_A yiQwi 

dt-\^jy^) 2/(r)e(r) /(r)|e(r)|; /(r) 

de 

dt " /(r) 

du 
dt 



If we take the energy integral H to have constant value h then the energy 
relation gives 

«gn(0W)(^)+7=0|:y (17) 

Similarly if we take the angular momentum pg to have constant value c, then 
the angular momentum relation takes the form 

c , , 

u 



Moreover, since prPe = uvf{r)\Q{r)\ the integrals Ii and I2 can be written as 

h = sm{f3e)\e{r)\f{r)uv - ^cos(/?^)e(r)|e(r)|/(r)V - ^ cos(/3^), 

^ (19) 

/2 = -cos(/?^)|e(r)|/(r)Mf -/3sin(/?^)e(r)|e(r)|/(r)V- ^sin(/?^) 
The system (1161) is no longer Hamiltonian, but (1171) defines a codimension one 



10 



invariant set 



yi{h) = {(r,9,u,v) e M^l r G [rN,rs] and the energy relation holds}. (20) 
which we continue to call energy manifold. 

System f|T6l) determines a vector field on M(/i) which is undefined when r = 
OT r = rs- We now consider the singularity at r = r^v. We will come back to 
the singularity at r = in Section [91 Let 

N = {{r,9,u,v) eM{h)\ r = Vn} 



which is the manifolds of states corresponding to collision. From the definition 
of M(/i), we see that 



N = {{r,9,u,v) e 



tn and 



2 I 2 

u +v 



27} 



and hence N is independent of h. Since 9 is considered modulo 27r, N is a two 
dimensional torus. 

Note that the vector field given by (ITB]) is not defined on N. The orbits ap- 
proaching N in a finite time are the collision orbits. However, we can scale the 
vector field in such a way that the new vector field can be extended to N. We 
accomplish this scaling by introducing a new time parameter r given by 



dr 



16(01 



dt. 



System ([5]) then, using Proposition [21 becomes 



r = vf{r) 
V = u^f'{r) 
9 = u 



2f{r)0{r) /(r)|e(r)| 



f22) 



u 



-uv 



{f'i 



r) + 



2/(r)e(r)y 

where the dot indicates differentiation with respect to r. For this new vector 
field, N is an invariant set. We call N the collision manifold. The flow on N is 
fictitious (i.e. it has no physical meaning), however, due to the continuity of 
solutions with respect to initial conditions, it gives informations about collision 
and near-collision solutions. 



7 Flow on the Collision Manifold 

Let k = y/KL. Since /'(t'at) = VKL > where \fKL = \(3\ by Proposition 
[21 Moreover, by the same proposition, /(r)9(r) ~k/P^ as r — and, near 
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the north pole, we have that 0(r) < 0. It is easy to verify that in the variables 
(r, V, 9, u) the points (r^r, ±^27, 6*, 0) are equilibria for system (l22l) . Therefore 
the sets 

S± = {{r,v,e,u) G N|t; = iy^} = {{r,v,e,u) G M^|(r, t;, ^, m) = (r^,±y2^,^, 

(23) 

are circles of rest points on the collision manifold N. At these points the 
linearized system has the matrix 



±fcv^ 

* ±P'^^/k 







1 







tV2^(A; - /?V(2A;)) 



(24) 



where * denotes an entry of the matrix that is not used in the computation 
of the eigenvalues. The corresponding eigenvalues are ±/c-\/27, ±/3^a/27//c, 
and T^{k- P^/{2k)). 

Using that f{r)Q{r) —k/jS"^ as r — and that 0(r) < near the north 
pole we see that the restriction of system (1221) to N yields the system 



V ■■ 
6 - 
u 



2k 



U 



U 



-UV 



2k 



where we have used the energy relation to simplify v. Introducing the angular 
variable x via 



M = Y 27 cos X 
V = 1/27 sin X 



(25) 



(26) 



one finds 

^=(k-^)=\^ 
de V 2A;y 2 ■ 

The solutions of this vector field are sketched in figure [2l 
In the resulting phase portrait w > so all orbits travel from the lower circle 
of rest points S~ to the upper circle S"*". Moreover, in the (6*, x) variables, the 
unstable manifolds of points on S~ are just straight lines of slope |/3|/2. The 
unstable manifolds at 9 = 6*, v = —y/2rf in N do not necessarily join up with 
the stable manifold at = f = ^/2rf. Only when we have 



dx 
d9 



k- 



0^ 
2k 



1 

2^ 



n 



1,2,3, 



(27) 
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-7T 



/2) 



Fig. 2. The flow on the collision manifold N of the generalized Kepler problem for 

we have this property. That is if = 1/n for a positive integer n, then each 
branch of the unstable manifold makes n circuits of N before rejoining the 
upper circle of rest points at the same ^- value. 
When 

2 



d0 [ 2k. 



n = 0,1,2,3, ... 



(28) 



2n + 1 

or equivalently when = 2/(2n + 1), the unstable manifold leaving 9 — 9* 
join up with the stable manifold at ^ = ^* + tt after making n + 1/2 circuits. 
In all other cases, the two branches of the unstable manifold reach distinct 
equilibrium points. 

Physical Interpretation. We now give a physical interpretation of the solution 
described above. As we mentioned earlier, the orbits on the collision manifold 
have no physical meaning, but give information about collision and near col- 
lision orbits. A near collision orbit makes n revolutions about r = before 
exiting at an angle which depends on j3. In the two special cases considered 
above the orbit either exits in the direction in which approached collision 
(1/3 1 = 1/n) or else in exactly the opposite direction (|/3| = 2/(2n + 1)). 
If 13 is not of one of these two forms, then nearby initial conditions will lead to 
quite different behavior near collision (this is the basic idea behind Easton's 
notion of topological regularization). In this case we cannot join orbits coming 
to collision with orbits leaving collision in a meaningful way so to make the 
resulting flow continuous. 



8 An Isolating Block About Collisions 

Let M be a smooth manifold and let : M x M ^ M be a flow on M. A 

subset N C M is called invariant if ■?/'(N,]R) = N. 

Definition 2. A compact invariant set N C M is called isolated if there exist 
an open set U containing N such that '?/'(a;, M) e U implies a; e N. The set U 
is called isolating neighborhood for N. 
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Now let B be a compact subset of M with non-empty interior and suppose 
that b = is a smooth submanifold of M. Denote by b"*" the set of ingress 
points of b, 

b+ = {a; G b| tpix, (0, -e)) n B = for some e > 0)}, 
by b~ the set of egress points of b, 

= {x e b| iIj{x, (0, e)) n B = for some e > 0)}, 

and by t the set of tangency points of b 

t = {a; G b| ip{x, 0) is tangent to b at x}. 

In general b"'' , b~ and t might be variously related but their union must always 
be b: points of b"*" leave B going backwards, those in b~ leave b going forwards, 
and the remainder must be in t. 

Definition 3. B is called an isolating block if t = b+ fl b~ (and t is a smooth 
submanifold of b with codimension 1). 

Definition 4. Let N be an isolated invariant set, and let B be an isolating 
block. Then B is said to isolate N if int(B) is an isolating neighborhood for 
N 

The following theorem was proved by Conley and Easton [1] 
Theorem 1. If N is an isolated invariant set, then there exists an isolating 
block which isolates N. If B is an isolating block, then there exists an isolated 
invariant set (possibly empty) which is isolated by B. 

We now want to define B in terms of a real valued function on M. Let / : 
M ^ M be a smooth map. We write 

r{x,t) :=/(V^(x,t)), 

and define 

j(x) = r (x, 0) and /(x) = /* (x, 0) 

where /* and /* denote derivatives with respect to time. The following lemma 
is proved by Wilson and Yorke \l2l (the symbol "D" denotes derivative). 
Lemma 1. Let J : M ^ [0,cx3), and let 6o > 0. Suppose that DI{x) ^ 
whenever < /(x) < 6o- Suppose also that /(x) > whenever < /(x) < 6o 
and /(x) = 0. Then N = /~^(0) is an isolated invariant set and /~^([0,5]) is 
an isolating block for N for each 6 G (0, ^o]- 
We now define the subsets of b that are asymptotic to N: 

a+ = {x G b+|V^(x, [0, oo)) C B}, 
a" = {x G b~|V^(x, (-00,0]) C B}. 

By definition b+ — a+ denotes the set of point in b with the property that 
the corresponding solutions enter in B without staying permanently there, i.e. 
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there exists a t > such that ipi^j t) ^ B. Thus we can define the time spent 
in the block for a point x G b"*" — a"*" by 

T{x) =inf{t > 0\ipix,t) ^B}. 

Note that ip{x, [0,T(x)]) G B and that iIj{x,T{x)) G b~. Now we define the 
map across the block 

\& : b+ - a+ ^ b" - a" : X tp{x, T(x)). 

Theorem 2 (Conley and Easton [1]). If B is an isolating block, then : 
b+ — a+ ^ b~ — a~ is a diffeomorphism 

Definition 5. An isolating block B is said to be trivializable if \Ef extends 

uniquely to a diffeomorphism from b"*" to b~. 

Trivializability is actually a property of an isolated invariant set: 

Lemma 2 (Conley[3])- Suppose that N is an isolated invariant set and that 

Bi and B2 isolate N. Then Bi is trivializable if and only if B2 is trivializable. 

Definition 6. Let B isolate N. Then N is said to be trivializable if B is 

trivializable. 

We now return to the generalized Kepler problems on spherical surfaces of rev- 
olutions. We take ip to be the fiow on the manifold 'WL{h) defined by equations 
( I2UI) . The invariant set N is given by fl?I]) . Define 

I -.Mih) ^R: {r,v,9,u) 



|e(r)| 

B{h,6) = {xe M{h) : /(x) < 6} 

Lemma 3. Given any h, there exists a 5o > such that B(/i, 6) is an isolating 
block for N whenever < 6 < Sq. 

Proof. The tangent space to M(/i) at the point x = (r, v, 6, u) is given by 

r h .... 

provided f , 0, u) 7^ (0, 0, 0, 0). Let x = (f, v, 9, u). Since DI(x) ■ x = 

j(r)2e(^)2 , we have that DI{x) 7^ if (m, v) 7^ (0, 0). The energy relation implies 
that 

|sgn(e(r))(u2 + v^) + 27I = 2\h\5 

where I{5) = 5 > 0. Therefore if 5q is small enough (i.e. 5q < 7/|/i|), then 
(■u,f) 7^ (0,0) whenever 5 < 5q. Hence DI{x) 7^ when < I{x) < 5q. Now 
using equation ( 122|) we see that 

sgn(9(r)) 1 



/(r)e(r)2 f{r)Q{ry 



15 



where the last equality follows from the fact that near the north pole G(r) < 0. 
Moreover (again assuming that 0(r) < 0, since we are near the north pole) 



/(r)0(r) 



5/2 



J(r)e(r; 



+ fir) + 



7 



/(r)e(r) 



If I{x) = 6 and if I{x) = 0, then \0{r)\ ^ = 6 and f = 0. Using the energy 
relation and Proposition [2] yields, 



/(r)e(r)2 
2 



/(r)e(r)2/'(r 



r(r)M + 7(^/'(r) + 
(M + 7)f(r)'- 



1 



2/(r)e(r) 



7- 



Since 0<5<5o<7/|^|"we have that hS + 'j > and the expression in 

square brakets is a convex parabola in f'{r). Moreover, since f'{rjq) = 

it follows that, if 6o is sufficientely small, then {h6 + 'y)f'{r)'^ — > 

and f'{r) > 0. Consequently J > whenever < 6 < 6o. Finally, note that 
N = /~^(0). Hence, by Lemma [T], B(/i, 6) is an isolating block for N. □ 

We now exhibit for the block B(/i, 5) the various subsets defined above. We 
fix h and choose < 5 < ^o, where 6o is given in Lemma O 

h = {xe M{h)\ e(r) = 6} 
h+ = {xeh\v <0} 
h- = {xeh\v>0} 
t = {x eh\ V = 0} 
a+ = {x e b+l M = 0} 
= {x E \ u = 0} 

Theorem 3. If the set N is a trivializable isolated invariant set for equations 
then \f3\ = 2/m, where m is a positive integer. 



Proof. Using the definition of B, we write 

1 



{{r,v,e,u) G 



|0(r)| 



6 and + = + 2h6}. 



The general form of the map across the block is 

: — — > — : (r, v, 9, u) — > (f, v, 9, u) 



we will use the first integral of the system to simplify this expression. Sup- 
pose {r,v,9,u) G b+ — a+ and let 'if{r,v,9,u) = {f,v,9,u). Recall that 
the energy h{r,v,9,u), the angular momentum c{r,v,9,u), Ii{r,v,9,u) and 
/2(r, V, 9, u) are first integrals of the problem. Then /i(r, f , ^, u) = h{f, v, 9, u), 
/i(r, V, 9, u) = /i(f, V, 9, u), /2(r, v, 9, u) = h^r, v, 9, u) and c(r, t>, 9, u) = c(f, v, 9, u). 
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Since r is constant in b it follows that f = r and the integral of the angular 
momentum yields u = u. From the energy integral and the definition of b"*" 
and b~ it follows that v = —v. Therefore we can write the map across the 
block as 

^ : b+ - a+ ^ b" - a" : (r, v, 6, u) (r, -v, u),u) (29) 

where \l/6i(^, u) is the third component of Here we are using {6, u) as coordi- 
nates on b+, so l/|G(r)| = 5 and v = —(2'j + 2h5~u'^Y^'^. Since the equations 
fl22l) are independent of 9 

^eiO,u) = e + T{u). 

The function F is defined for all u such that < u"^ < 2'~f + 2h6. By symmetry, 

T{-u) = -T{u). (30) 

Assume that B is trivializable. Then \1/ extends to a continuous map b+ — > b^. 
Thus 

9 + F(0+) = 9 + F(O-) + 27im. 
where m is an integer. By (150]) . F(0— ) = — F(0+). Hence we must have 

F(0+) = irm. 

The number F(0+) can be computed using geometric methods. Consider a 
point xo G a"*". The orbit through xq is the stable manifold of a point s~ in 
S~. Now let a; G b"*" be close to xq. the orbit through x follow closely the stable 
manifold of s~. 

In Section 7 we studied the flow on N where we found out that the unstable 
manifolds of points on S~ are just stright lines with slope |/3|/2. We are inter- 
ested in the branch of the unstable manifold of s~ for which u > 0. Therefore 
we take — 7r/2 < x < ^/2- Write 

s- = C5,9*,0,-^), 

where 1/|G(5) = 6\. Then the unstable manifold of s~ is exactly the stable 
manifold of the point 

s+ = (5,r + 27r/|/3|,0,V^). 

We now can determine F(0+). The orbit through x first follows the stable 
manifold of s~, then follows the unstable manifold of s~, which coincides with 
the stable manifold of s"*", and finally follows the unstable manifold of s+. 
Therefore, the change in 9 along the orbit approaches the difference 

in 9 between s~ and s"*". This difference is 27r/|/3|. Hence F(O-I-) = 27r/\(3\. 
Since we also have that F(0+) = irm we obtain = 2/m. Since \(3\ > 0, m 
is positive. This completes the proof. □ 
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We now want to prove that N is a trivializable isolated invariant set for equa- 
tions (122|) . In order to do that, unhke in the work of McGehee we will use 
the additional integrals of motion of the system. 

Theorem 4. If = 2/m, where m is a positive integer, then the set N is a 
trivializable isolated invariant set for equations fl22|) . 

Proof. Recall that 

b = {(r,v,9,u) e M^l— ^ = S and u'^ + v'^ = 27 + 2M}. 

|e(r)| 

Suppose (r, v, 9, m) G b"*" — a"*" and let \l/(r, v, 9, u) = (f, v, 9, u). Then, in order 
to show that N is trivializable, we must prove that the map across the block 

: b+ - a+ ^ b" - a" : (r, v, 9, u) (f, v, 9, u) (31) 

extends to a diffeomorphism of b+ onto b^. Here, as shown in the proof of 
Theorem 3, m) is the third component of and we are using {9,u) as 
coordinates on b"*", so l/|G(r)| = 5 and v = —(27 + 2hS — -u^)^/^. 
Since Ji and I2 are integrals of motions we have Ii{r,v,9,u) = Ii{f,v,9,u) 
and, l2{'r, v, 9, u) = /2(^, v, 9, u). Therefore we have 

Ji{r,v,9,u) = h cos (39+l2sm [39 = I cos{(3{9-9o)) = -/90(r)|e(r)|/(r) V-^ 
and 

J2{r, V, 9, u) = h sin (39 - h cos f39 = I sm{(3{9 - 9o)) = \e{r)\f{r)uv, 
where 

Ii = I cos p9o, l2 = I sin p9o (32) 

and / = \J Ii + li- Consequently Ji(r, v, 9, u) = Ji(f, v, 9, u) and J2{'r, v, 9, u) = 
—■hir, V, 9, u). This gives the system of equations 

I cos{p{9 - 9q)) = I cos{(3{9 - 9o)) 
I sm{(3{9 - 9o)) = -I sm{(3{9 - 9o)) 

where 9q = 9f){f, v, 9, u). Since 6*0 = ^0 the general solution of the system above 
is 

m, = e = -9 + 29^ + ^. (33) 

We now want to write '^q{9,u) in the form '^q{9.,u) = 9 + r(M) used in the 
proof of Theorem 3. First we can rewrite equations f|T9|) as 

/i = /sin[/3(0-C)], l2 = -Icos[[3{9-C)l (34) 

where cos(/3C) = (|e(r)|/(r)Mf )// and sin(/?C) = (/3e(r)|e(r)|/(r) V+7//?)//. 
If we suppose {r,v,9,u) G b+ — a+ then l/|9(r)| = 6 and if we use {9,u) as 
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coordinates on b+ we have 

-fclM(27 + 2M-M2)l/2 



cos(/5C) 
sin(/?C) 



(A;2(27 + 2h6 - u^)u^ + {k^u^ + ^)2)i/2 



2 _i_ 7^ ("^^) 



(A;2(27 + 2M - m2)„2 + (^^m^ + ^)2)i/2^ 



where k\ and /c2 are constants such that fci = |0(T')|/(r) and ^2 = /30(r)|6(r)|/(r)2M2+ 
7//5, when l/|e(r)| = 5. 

Note that the right hand side of the second equation of system fl35|) is always 
positive. It follows that < /3C < tt. Therefore the function C, is defined by 

, 1 / -A;im(27 + 2M-m2)V2 

= — arccos ' 



/3 \ (A;? (27 + 2M - m2)m2 + (^2^2 + |)2)i/2 

and it is of class C°° since composition of functions (the inverse cosine 
function is of class C°°). 

From equations fl32|) and flMl) we obtain the following system of trigonometric 
equations 

cos/^^o = sin[/3(0 - C)], sin/?^o = -cos[/?(^ - ()]• 

The general solution of the system above is 6*0 = 6* — C — ^ + We choose 
n = and we define as 

° ^ 2/5 

Substituting the above expression for 6*0 in fl5^ we obtain 

vi/, = ^ + r(..) = ^-2C ^^"^ 



where r(-u) = —2^ ~ f + ^^"-^ ^ integer to be determined. 
We now need to determine the value of k imposing continuity at m = 0. Recall 
that r(0+) = vrm = ^. Since /?C(0+) = f we have 1(0+) = ^(fc - 1) = ^. 



/3 

Consequently, ifu>OA;=l + ||. Since r(O-) = -1(0+) and r(O-) = 
^(/c — 1) = Pi we find that, if m < 0, A; = 1 — -p. Therefore we can define 
r(ti) as 



-2C-i + f(l + ^) ifn>0. 



We can now show that the map across the block ^ extends to a diffeomorphism 
^ from b+ to b~. Let 

^ : b+ ^ : (r, v, 0, u) (r, -v, 9 + T{u), u), 
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where 

l-2C-| + t(l + M) if^>0 
be the extended map. Such map is continuous by construction. Moreover, since 
({u) is differentiable, is a differentiable map ^. Moreover, let {f,v,6,u) G 
b~, then : — b"*" is defined as follows 

§"^(f, V, 9, u) = {f,9 - T{u),u, -v) 

is the inverse of the map The differentiability of \E'~^ follow immediately. 
Therefore extends to a diffeomorphism, this concludes the proof. □ 

9 Block Regularization of the Vector Field 

In the previous section we worked with system (122|) . for which N is an invariant 
set. We now turn to our original problem: determining if equations ([5]) can 
be regularized. We consider the equivalent system (fT6ll . The set N given by 
( l2Ti) is the set of singularities of equations ( fT6l) . i.e. it is the set where the 
vector field fails to be defined. On 'WL{h) — N the orbits for the two sysyems 
are identical: only the parametrization is different. We now turn, following 
Easton, to the definition of block regularization. 

Let M be a smooth manifold, let N be a compact subset of M, and let F be 
a vector field on M — N, where, in this section, N is the set of singularities of 
the vector field F. Let be the fiow on M — N given by F (we do not require 

to be defined for all t). 
Let B be a compact subset of M with non-empty interior, and suppose that 
b = dB is a smooth manifold which does not intersect N. Let b+,b~ and t 
be defined as in the previous section. Let the definition of isolating block be 
also the same. Let 0{x) denote the orbit through x, namely 

0{x) = t) is defined}. 

Definition 7. An isolating block B is said to isolate a singularity set N if 
N C int(B) and if 0{x) ^ B for all x G B - N. 

The subsets a"*" and a~ are the same as before, except that now we must allow 
for solutions which are not defined for all t. Therefore 

a^ = {x G h'^\(j){x,t) G B for alH > for which 0(a;,t) is defined}, 
a^ = {x G b^|0(x, t) G B for alH < for which 0(x, t) is defined}. 

The map $ : b+ — a+ ^ b~ is defined in exactly the same way as the map \E' 
and we have that $ : b+ — a+ — > b~ — a^ is a diffeomorphism. We also have 
the same definition of a trivializable block B. 

Definition 8. The singularity set N is said to be block regularizable if there 
exists a trivializable block B which isolates N. 
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7 



Fig. 3. Cone as quotient of a plane sector 

We now return to our original problem. Note that whether a certain set is or is 
not an isolating block is independent of the parametrization of the flow. The 
map across the block is also independent of the parametrization. Therefore, 
B{h, 6) is an isolating block for system ([T6|) if and only if it is an isolating 
block for system fl22l) . and $ = Hence B(/i, 6) is trivializable for f|T6l) if and 
only if it is trivializable for system fl22|) . 

Theorem 5. The singularity set N (i.e. the singularity at the north pole) 
for system flTB]) is block regularizable if and only if \j3\ = 2/m, where m is a 
positive integer. 

We can now quickly discuss the singularity set at the south pole for system 
(fT6l) . In this case the situation is much simpler and the singularity is always 
block regularizable. Let N5 be such singularity set, i.e. 

N5 = {{r,9,u,v) e MV = rs} 

then one can prove the following 

Theorem 6. The singularity set N5 for system ( |T6l) is block regularizable. 

Proof. We give a sketch of the proof. One can construct an isolating block 
repeating the proof of Lemma [3] with some minor changes. In this case the sets 
a"*" and a~ are empty. This is a consequence of the fact that, as it can be seen 
from the energy relation ([T7]) . no solution ever reaches the singularity set N^. It 
follows from Theorem [2] that the map across the block \1/ : b+— a+ b^ — a~ is 
a diffeomorphism. Moreover, since a"*" and a~ are empty, \I' trivially extends to 
a diffeomorphism from b"*" to b~ . Therefore the isolating block is trivializable. 
Moreover, as we noticed above, a region is an isolating block for system flTBl) 
if and only if it is an isolating block for fl221) . and $ = \I'. Since the isolating 
block can be constructed so that it isolates N5 it follows that N5 is block 
regularizable. □ 

10 A Note About Orbifolds of Revolutions 

Spherical orbifolds of revolutions are easily identifiable by their tangent cones 
at the poles. Namely, the tangent cone at the pole must be isometric to the 
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metric quotient of the flat plane by a finite cyclic group of rotations fixing 
the origin. Note that the tangent cone at the pole is generated by rotating 
the tangent line to the profile curve at the pole about the axis of rotation. If 
the cyclic groups at the poles are of different orders, the orbifold is commonly 
referred as bad since it will not arise as a quotient of a Riemannian 5^ by a 
finite cyclic group of isometries [TTj . 

In general a fiat right circular cone with vertex angle is obtained by iden- 
tifying the edges of a plane circular sector of angle 7. The relation between 
(p and 7 is easily computed: 7 = 27r sin 0. See figure O Thus if the tangent 
cone at a pole of a spherical orbifold of revolution is isometric to M^/Z^, then 
7 = 27r/n for a positive integer n. 
We are now in a position to prove the following. 

Theorem 7. The singularity set N for system (fT6|) is block regularizable if 
the surface of revolution is an orbifold of revolution with constant curvature. 

Proof. Let (j)N, resp. (ps, be the angle between the profile curve 7(r) = {f{r),0,g{ 
and the axis of rotation at r = r^, resp. r = rs- For an orbifold of revolution 
we must have sin (p^ = 1/n and sin (ps = ^/k for some positive integers n and 
k. 

On the other hand siiKpN = fir^) = \(3\. By Theorem 5 the singularity N 
for system f|T6|) is block regularizable if and only if \(3\ = 2/m, where m is a 
positive integer. 

Therefore for an orbifold of revolution we have m = 2n. A similar argument 
also shows m = 2k. This proves that the singularity N is always block regu- 
larizable on orbifolds of revolution with constant Gaussian curvature. □ 
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